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Transistors constitute the backbone of modern day electronics. Since their advent, researchers have been 
seeking ways to make smaller and more efficient transistors. Here, we demonstrate a sustained amplification 
of magnetic vortex core gyration in coupled two and three vortices by controlling their relative core 
polarities. This amplification is mediated by a cascade of antivortex solitons travelling through the dynamic 
stray field. We further demonstrated that the amplification can be controlled by switching the polarity of the 
middle vortex in a three vortex sequence and the gain can be controlled by the input signal amplitude. An 
attempt to show fan-out operation yielded gain for one of the symmetrically placed branches which can be 
reversed by switching the core polarity of all the vortices in the network. The above observations promote the 
magnetic vortices as suitable candidates to work as stable bipolar junction transistors (BJT). 
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There has been a revolution in the study of inhomogeneous and non-trivial magnetic nano-structures such as 
magnetic vortices and antivortices due to their suggested applications in magnetic data storage, magnetic 
random access memory 1 " 4 , magnetic logic 5 and information processing devices 5 . The presence of rotational 
asymmetry makes them very attractive candidates for studying the interaction between the local magnetization 
and externally applied magnetic fields or spin polarized currents 6 " 16 . Magnetization can align itself parallel to the 
edges of a thin ferromagnetic nano-disk in order to reduce the surface charges. This gives rise to the clockwise 
(CW) or counter-clockwise (CCW) chirality of the vortex. The core of this vortex, which is a few nanometres in 
diameter, gets pushed out (p = 1) or into (p = — 1) the nano-disk's plane as the exchange interaction favours a 
parallel arrangement. This is defined as the polarity of the vortex. 

Magnetic vortices can be brought to gyration by the application of magnetic fields or spin polarized cur- 
rents 116 " 20 . In addition to the external forces, the moving vortex core experiences another internal force arising 
from the demagnetizing field of the non-equilibrium magnetization pattern. This force is perpendicular to the 
vortex core's velocity and leads to its gyrotropic motion. For large amplitude excitation, vortex core switching 
occurs along with the creation and annihilation of new vortex and antivortex pair 21 ' 22 . For small amplitude 
excitation, vortex core motion remains in the linear regime 23 . 

The magnetization dynamics resulting in the gyrotropic vortex core motion is described by the Landau- 
Lifshitz- Gilbert equation 24 ' 25 . In the linear regime, the vortex core equation of motion can be derived from the 
Thiele's equation 23 . The CW or CCW sense of vortex core gyration direction is solely controlled by its polarity. In 
the linear regime, vortex core can be described by a harmonic oscillator model 26 . Consequently, magnetostatically 
coupled vortex gyration can be considered as coupled oscillators. Therefore, one expects mutual energy transfer 
and a consistent phase relation between the gyrating vortices 27 ' 28 . Logic operations based on magnetic vortex state 
networks have been demonstrated experimentally via the vortex gyration mediated information signal transfer 
mechanism 29 . 

Subsequently, the mutual transfer of energy between magnetostatically coupled vortices, where one of the 
vortices is locally excited, is extremely important for microwave communication and logic operations. In this 
regard, the parameters like the signal transport rate and efficiency are the key factors in determining the device 
performance. Vortex gyration transfer rate and energy attenuation coefficients have been calculated by analytical 
method and micromagnetic simulations 30 . Stimulated vortex gyration based energy transfer between spatially 
separated dipolar coupled magnetic disks has been observed by time resolved soft x-ray microscopy 13 . The rate of 
energy transfer is found to be determined by the frequency splitting caused by the dipolar interaction between the 
vortices 29 . 

This energy transfer efficiency may depend on several factors such as the frequency of the exciting field pulse as 
compared to the gyration frequency of the vortex core, the amplitude and nature of the exciting pulse, the distance 
between the vortices and their relative polarity. Until recently 13 , the energy transfer efficiency was found to be well 
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below 100% 27 ' 28 and no amplification has been reported thus far. It 
has also been reported that the interaction strength between coupled 
vortices is maximum when their core polarities are opposite 13 , 
although, higher interaction strength is not the sufficient condition 
for higher transfer efficiency. For higher amplitude input, the vortex 
motion enters the non-linear regime and vortex core switching 
occurs; making it impossible to increase and maintain a large ampli- 
tude output and a constant phase relation between gyrotropic 
motion of both the vortices. On the other hand, if the input signal 
is very weak and the frequency is close to the gyration frequency of 
the vortex core, the amplitude of the response gradually increases 
indicating that the core switching may occur at some point, which is 
not desirable for device application. 

In our results we show that off- resonant signals 31 of lower ampli- 
tude can be used to design suitable transducers with isolated vortices, 
which will be required to convert other kinds of signals (e.g. a rotating 
field) to vortex core gyration. In the case of a pair of magnetostatically 
coupled vortices, if a signal is applied to only of them then the other 
one shows a greater core gyration i.e., amplification when the core 
polarities are opposite. Antivortex solitons moving through the stray 
field are held responsible for this behaviour. We postulate some rules 
regarding their dynamics and use them to mimic transistor-like 
operations of switching and amplification with a chain of three vor- 
tices. Furthermore, we attempt to couple the output of this three 
vortex chain to two symmetrically placed daughter chains in an 
attempt to demonstrate a fan-out operation. However, the antivor- 
tices involved in the dynamics favoured one branch over the other 
resulting in a higher level of asymmetry - one of the branches 
received more power than the other. 

Results 

We use permalloy (Py: Ni 80 Fe 2 o) with negligible magneto-crystalline 
anisotropy in the form of a 40 nm thin disk of diameter 2R = 200 nm 
to ensure a stable vortex structure 32 . The darker shade in Fig. 1 (a) 
represents such an isolated vortex. Figure 1 (b) shows a pair of 
coupled vortices whose centre to centre distance is a = 250 nm. A 
chain of three vortices, with the same centre to centre distance a, has 
also been studied with different orientations of polarity. Spatially 
averaged X-component of magnetization (m x ) has been used as an 
indicator of core displacement away from their equilibrium posi- 
tions. The square of the amplitude of Fourier transform of (m x ) (t) 
(with respect to t) - also known as the energy spectral density (ESD) - 
shows the peaks in vortex core dynamics as a function of frequency. 

Isolated magnetic vortex. Figure 1 (c) shows a plot of (m x ) vs. time 
and Fig. 1 (d) shows the associated ESD (in decibel) for the single 
vortex excited by a broadband signal. The gyrotropic mode is 
observed at frequency/ = f 0 ~ 1.27 GHz. Higher frequency modes 
associated with the generation of spin-waves 33 are also observed. 
Here, we concern ourselves with frequencies / < f 0 while using 
signals that are rotating CW or CCW in the plane of the magnetic 
vortices. With up polarity p = 1, CCW signals are known to produce 
greater gyration 10 which can lead to polarity switching. This can be 
useful in terms of data storage 31 . However, if the polarity switching is 
somehow avoided, one can use this to create a suitable transducer for 
appropriately rotating signals. To that end we can use signals with 
lower amplitude. Signals with off- resonant frequency/ </ 0 , should 
be used to reduce the convergence time. As shown in the Supple- 
mentary Fig. Fl, anharmonicity of the gyration dynamics results in a 
beating frequency when off-resonant signals are used. A trade-off 
between convergence time and beating frequency needs to be further 
explored as a design consideration. 

Coupled magnetic vortices pair. We have examined the transfer of 
energy from one vortex to another in terms of their core gyration 
amplitude (measured in terms of (m x )) when excitation is only given 




Figure 1 | Isolated vortex and coupled vortex pair. Dark regions represent 
the 40 nm thick (a) isolated and (b) coupled magnetic vortices each of 
diameter 2R = 200 nm. The centre to centre separation in the case of 
coupled vortices is set to a = 250 nm. (c) Time evolution and (d) 
corresponding energy spectral density of (ra x ) in response to the signal H%. 

to one of them. Here, (m x ) is computed for both the vortices sepa- 
rately. The dynamics was examined with all sixteen combinations of 
polarity and chirality of the two vortices. When a small external bias 
field is applied along the X-axis, the vortex cores may move up or 
down along the 7-axis. This changes their separation and causes 
magnetic surface charges to appear on the vortex boundaries; 
consequently affecting the strength of their magnetostatic 
coupling 34 ' 35 . Hence, in the presence of a bias field, if both the 
vortices have the same chirality, their coupling will remain 
relatively unaffected, than when they have different chiralities. This 
phenomena can be used to affect a chirality dependent dynamics and 
signal transmission. However, in the absence of an external bias, we 
observed that chirality does not play any role towards enhancing the 
asymmetry in dynamics, as described below. Thus, henceforth all 
vortices in this work have CCW chirality. Furthermore, observable 
changes only appeared to occur between cases with similar and 
opposite polarities 14 . Mediated by several factors 8 ' 36 " 39 , the resonant 
frequencies of a pair of vortices can differ from that of an isolated 
vortex. However, in this study, we used an excitation signal rotating 
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Figure 2 | Gyration amplification. ESDs of left and right magnetic vortices 
shown in the insets with ((a) and (c)) similar and ((b) and (d)) opposite 
polarities, (a) and (b) show the results for a signal amplitude of 0.5 mT 
while (c) and (d) show those for an amplitude of 1.5 mT, both rotating 
CCW at/=/ 0 . 

at frequency f 0 , which is applied only on the left disk (Fig. 1 (b)). 
Figures 2 (a) and (c) show the results for the case when both polarities 
are up (pip 2 = 1) and Figs. 2 (b) and (d) show those when left core is 
up and right core is down (pip 2 = — 1). Figures 2 (a) and (b) 
correspond to a signal amplitude of 0.5 mT and Figs. 2 (c) and (d) 
correspond to a signal amplitude of 1.5 mT. 

With increase in the signal amplitude from 0.5 mT to 1.5 mT, we 
see a splitting 40 of the peak. Furthermore, when pip 2 = — 1, more 
energy is transmitted and stored in the right vortex. For a signal 
amplitude of 0.5 mT and p x p 2 = 1, the left vortex, exhibits 
60.72 dB of ESD at its gyro tropic mode while the right vortex exhibits 
68.43 dB. While for p x p 2 = -1, these values become 30.99 dB and 
53.69 dB, respectively. Thus, the difference in ESDs of gyrotropic 
modes of left and right vortices increases by about 15 dB. When 
the signal amplitude is increased to 1.5 mT, these values become 

81.66 dB (left vortex) and 79.76 dB (right vortex) for p x p 2 = 1 and 

46.67 dB (left vortex) and 72.45 dB (right vortex) for p x p 2 = -1. 
Here (for p x p 2 = —1), the ESD (at / = f 0 ) of the right vortex is 
25.78 dB greater than that of the left vortex. This shows that opposite 
core polarity facilitates amplification of signal transfer, which is fur- 
ther enhanced by the input signal amplitude. An increase in about 
9.5 dB in signal power (from 0.5 mT to 1.5 mT) has caused the 
difference in gain to increase by 12.69 dB. The dependence of this 
relative amplification on the strength of the input signal is essential to 
mimic the transistor operation where the base current controls the 
amplification. 

These observations, along with the ones made for an isolated 
vortex, testify to the existence of anharmonic and asymmetric 
dynamics present in the vortex core gyration, which cannot be 
explained by solutions of the Thiele's equation with linear approx- 
imations 23 ' 41 ' 42 ; even if vortex core deformation 43 is taken into 
account. Although, the type of amplification described here has 
not been seen before, asymmetry based on polarity in terms of energy 
transfer rate and efficiency has been observed experimentally. 
Stronger or weaker stray field coupling can affect the rate of energy 
transmission 13 , but it does not guarantee the observed asymmetry in 
general and amplification in particular. 




Figure 3 | Stray field distribution showing the path of travelling 
antivortex packets for p\p 2 — 1 (a) and p x p 2 = —1 (b). With time the 
packets shift their path from the dashed to the solid lines. 

One may draw an analogy of this observation with a driven double 
pendulum made of identical pendulums. In this case, when the 
driving frequency is same as the eigenfrequency of the isolated pen- 
dulum then the top pendulum mass becomes stationary while the 
bottom pendulum mass moves with an amplitude twice as that of the 
driving amplitude. This may be considered as an infinite amplifica- 
tion (although, the angular displacement of the top pendulum will 
still not be zero). However, for a coupled pendulum, where two 
pendulums are connected with a spring and dynamics is not pinned 
at any end (this system is closer to the coupled vortices presented 
here in terms of the underlying equations of their dynamics), no such 
amplification occurs when the driving frequency is the same as the 
eigenfrequency of the individual pendulum 44 (see Supplementary 
Note Nl). Either way, a direct comparison of the energy transfer 
mechanism of coupled magnetic vortices with that of coupled mech- 
anical oscillators is difficult due to the presence of additional para- 
meters in the former case. The gyrovector, which serves as the inertia 
of the vortex core 7 can switch direction with core polarity, resulting 
in the amplification observed here. Hypothetically, this will be com- 
parable to obtaining an amplification in one of the coupled pendu- 
lums by changing its inertia from J to —I. An analytical model is yet 
to be developed to describe this phenomenon; but, this is outside the 
scope of this paper (see Supplementary Note Nl). Here, we consid- 
ered a numerical approach and calculated the temporal evolution of 
the stray field and discovered that packets of antivortex structures 
travelling through the stray field mediate the transfer of energy 
between the two vortices. The path of these antivortices are shown 
in Fig. 3 for both polarity combinations when the excitation signal 
amplitude is at 1.5 mT. As time progresses the path shifts from 
dashed to solid lines. 

As seen from Fig. 3 (a), a single antivortex packet moves in a closed 
path for pip 2 = 1. This packet collides with other antivortex struc- 
tures which originate at the boundaries of the nanodisks. We under- 
stand that the antivortex packets discussed here are not particles in 
the true sense and their apparent 'collision' is only a result of the 
interacting stray fields. This interaction (or collision) is also shown in 
the Supplementary Movie Ml. It is only during this collision that 
these antivortex 'solitons' ('soliton' has been used loosely here to 
describe even those short-lived antivortex packets which do not 



SCIENTIFIC REPORTS | 4:4108 | DOI: 1 0.1 038/srep041 08 



3 



possess consistent form for significant duration) are allowed to 
change their size (local field distribution) significantly 45 . When it 
collides with the left disk, which initially has more gyrotropic energy, 
it becomes smaller and when it collides with the right disk, it becomes 
larger. This indicates that there is an inverse relation between the size 
of the antivortex and the gyrotropic energy that it can transfer. As 
time progresses, and the amplitudes of the gyrating vortex cores 
become comparable, the path of this bouncing soliton becomes smal- 
ler and it moves to a location shown by the solid line ® in Fig. 3 (a). 
The soliton itself does not change greatly in size after this point 
indicating no significant transfer of energy. 

Figure 3 (b) shows that more than one antivortex solitons are 
involved in the energy transfer for the case when p x p 2 = — 1. On a 
given vortex boundary, as one soliton gets terminated, another one is 
created. This creates a cascade of solitons, which vary in size (see 
Supplementary Movie M2). At first, the dashed lines mark the path of 
this cascade with the large arrowheads showing where a soliton gets 
localized. Branches are numbered from ® to ®, in an order such that 
the path of the new soliton is shown by the next branch. The path of 
this cascade changes gradually with time as well. Specifically, the 
number of rebounds between the boundaries of the vortices (the 
length of branch ®) may vary quite often. However, it is observed 
that as branch ® terminates, the remaining solitons, which follow 
branches © and © are relatively smaller. When branch © termi- 
nates, a soliton of the same size as the first one (which traversed 
branch ®) emerges from the right vortex to trace the final branch 
®. Thus the right vortex gains gyrotropic energy in the beginning of 
the dynamics. This cascade occurs twice every cycle, suggesting that 
signal transfer rate or efficiency can be controlled not only by the 
saturation magnetization 13 , but also by the frequency of the signal; 
and further optimization of signal transfer efficiency by tuning the 
dimensions of the coupled vortices is possible. As time progresses, 
the cascade starts to occur along the solid lines ® to ®, shown in 
Fig. 3 (b). When the gyration amplitude of the right vortex becomes a 
certain degree greater than that of the left one, we notice that the 
soliton, which was traversing the dashed branch ® earlier, now starts 
from the boundary of the left vortex. However, it is deflected back by 
another soliton, which emerges from the right vortex - much like an 
electron or hole charge carrier being prohibited from crossing the 
depletion layer of a junction diode. We can turn this amplification 
'on' or 'off simply by switching the polarity p 2 ; but it may be tech- 
nologically more desirable to have this control via a third vortex. 

Magnetic vortex transistor (MVT). In order to examine this 
transistor-like behaviour, we now add another vortex towards the 
right of the vortex pair shown in Fig. 1 (b) to form a three vortex 
sequence with polarities (from left to right) pi, p 2 andp 3 , which take 
values of 1 or — 1 denoting up or down polarities. In the previous sub- 
section, we identified relative polarity as the source of the observed 
amplification. Hence, here we study only the four cases withp x = 1 
(up), p 2 = ± 1 andp 3 = ±1. Chirality in all cases is CCW. Signal is 
applied to the left vortex only. The ESDs for these cases around 
frequency f 0 are shown in Fig. 4 (as shown in the insets, the 
excitation is given to shaded vortices only). 

Splitting can be observed in a few cases in Fig. 4. Unlike the 
splitting seen with increase in signal amplitude (see Fig. 2 (c)), which 
happens due to inherent non-linearities of the dynamics 40 , here it 
occurs for a different reason: an increase in the number of vortices 
leading to an increase in the number of permutations of couplings in 
the system 46 . Below, we consider any difference in ESD at the signal 
driving frequency off = f 0 only. 

As seen in Figs. 4 (a) and (c), the transmission efficiency is equi- 
valent for a persistent signal in both cases: (1,1,1) and (1,-1,1); with 
the latter faring slightly better. Although, a third vortex was added in 
the chain, a gain of 12.84 dB (between right and left most vortices) is 
observed in Fig. 4 (d). Also, transistor like switching is observed 
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Figure 4 | Magnetic vortex transistor. ESDs of left and right magnetic 
vortices with (pi,p 2 >p3) equalling (a) (1, 1, 1), (b) (1, 1, -1), (c) (1, -1, 1) 
and (d) (1, — 1, — 1) as shown in the respective insets. A 1.5 mT signal 
rotating CCW at frequency f 0 is applied only to the left (shaded) vortex, (e) 
Gain B, versus logarithm of signal amplitude h 0 . 

clearly with the three vortex sequence considered here when chan- 
ging from p 2 = — 1 (high base current) (Fig. 4 (d)) top 2 = 1 (low base 
current) (Fig. 4 (b)) changes the difference in signal levels of the right 
vortex (collector) from 12.84 dB to — 15.71 dB. We define the gain B 
in Eq. 1 as below: 



5^ESD 3 (f 0 )-ESD 1 (f 0 ) 



(1) 



Here ESD^o) and ESD 3 (f 0 ) are ESD at / = f 0 for left and right 
vortices, respectively. 

We further checked if this gain B, also depended upon the input 
signal amplitude h 0 . Figure 4 (e) shows a plot of B versus h 0 for h 0 = 
(le-6, le-5, le-4, le-3, 5e-4, 0.25, 0.5, 1, 1.5, 2) mT. Left vortex's core 
reversed for h 0 = 3 mT; and hence, we limit ourselves to 2 mT. For 
lower values of h 0 , the gain appears to be constant at B = B active ~ 
14.8 dB. This is reminiscent of a BJT operating under small-signal 
conditions 47 . At higher signal strength, the gain B, no longer remains 
constant. This indicates that like other electronic transistors, our 
'magnetic vortex transistor' is also susceptible to non-linear distor- 
tion. The maximum value of gain B = B max ~ 15.21 dB is seen for a 
signal strength of h 0 = /z max = 0.25 mT. 

We further investigate the roles of the stray field antivortex soli- 
tons on the transistor-like operations described above. We begin by 
analysing the temporal evolution of the stray field for cases where (pi, 
p2> Pi) equals (1, —1, 1) and (1, —1, —1). The same is shown in 
Supplementary Movies M3 and M4, respectively. Figures 5 (a) and 
(b) summarize the path of the solitons involved. Polarity dependent 
transient gyrotropic energy transfer between any two neighbouring 
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Figure 5 | Vortex chain. Stray field distribution in the cases where (pi,p2>p3) equals (a) (1, — 1, 1) and (b) (1, — 1, — 1). The path of antivortex packets after 
the dynamics has stabilized is marked with solid lines. 



vortices here too occur in the same manner as shown by the dashed 
lines in Fig. 3. These lines are omitted in Fig. 5 for clarity. Solid lines 
show the approximate paths that the solitons follow after the 
dynamics had become relatively stabilized. For (p l5 p 2 , p 3 ) = (1, 
— 1, 1), where no amplification (B < 0 dB) is observed, the cascade 
of solitons form a large oval loop around the central vortex. Energy 
appears to be transferred during collisions at ® and ©. This creates a 
closed feedback loop directly between the left and the right vortices. 
The solitons skip along the boundary of the central vortex on several 
occasions in order to aid their own cascade. Most importantly, we 
note here as well that solitons in the bottom half of the loop (right to 
left vortex) appear larger (lesser energy) than those traversing the top 
half (left to right vortex). Thus the flow of energy still occurs from the 
left vortex, which is excited externally, to the right one. However, an 
amplification (B > 0 dB) is observed for (pi,p 2 , p 3 ) = (1, — 1> — 1) 
and Fig. 5 (b) sheds some light on this crucial finding. Here, the path 
of antivortex solitons between the left and the central vortex does not 
change greatly from its early transient stage. Here too, energy is 
transferred at the sites ® and ©, essentially from the left vortex to 
the right one. However, unlike in Fig. 5 (a), no feedback loop, and 
thus no energy rebalancing is present here. This leads to a unidir- 
ectional flow of energy as determined by the cascade of antivortex 
solitons. The right vortex core, thus builds up gyrotropic energy until 
its drag (or dissipation) matches the power influx. Thus an amp- 
lification of the gyrotropic mode of the right vortex is observed in this 
case. One can simplify the dynamics for the two cases analysed above 
by considering the central vortex as an efficient medium and taking it 
out of the picture. Then we can see that amplification was observed 
when p!p 3 = — 1. However this amplification can be controlled by 
switching the polarity of the middle vortex (p 2 ), similar to what is 
done by switching the base current in a BJT. 

Fan-out. In an attempt to demonstrate a fan-out operation, which 
may support the development of more complex circuits, we placed 
two more MVTs symmetrically above and below the original MVT as 
shown in Fig. 6. Same material and structural parameters as before 
were used here. Power was given only to the left vortex of the original 
MVT. To our surprise, amplification was seen in only one of the 
branches. As shown in the inset of Fig. 6, the right vortex of upper 
branch in this network received about 15 dB more power than the 



lower branch. When all the core polarities in this network were 
reversed, the lower branch received the greater power by the same 
amount. 

The cause of this asymmetry is the fact that the solitons do not split 
during a fan-out. Also, the antivortex seen between the first two 
vortices in Fig. 5 (b) goes on directly to hit the upper chain as marked 
in Fig. 6. Thus, further study regarding the implementation of a fan- 
out is warranted by looking at the path of these antivortices in dif- 
ferent network configurations. 

Discussion 

We numerically examined the polarity dependent asymmetry and 
non-linearities in vortex dynamics. Cases presented in this paper 
included isolated vortices and coupled two and three vortex 
sequences. We particularly examined the dynamics for gain in the 
transfer of gyrotropic mode power from one vortex to another. To 
start with, we describe the design considerations in creating a trans- 
ducer which can convert power from a rotating field signal to gyro- 
tropic power. Best results were observed when the driving signal 
frequency was close to the eigenfrequency of the transducer. In the 
case of coupled pair of vortices, when an excitation signal is applied 
to only one of the vortices then considerably more energy (a max- 
imum gain of 27.68 dB) is transferred and stored in its neighbouring 
vortex if it has the opposite polarity. We further observed that this 
amplification of energy transfer can be extended over three vortices 
for a particular case of (p l5 p 2 , p 3 ) = (1, —1, —1). We interpreted 
these remarkable observations using the temporal evolution of stray 
magnetic field and observed that antivortex packets moving through 
the stray field were accountable for the observed amplifications. The 
rules, which we postulated based upon the motion of the antivortex 
packets (or 'solitons') can also successfully explain the previous 
experimental observations in greater detail. We hope that further 
study of these solitons will aid the research community in creating 
a better analytical model which can predict such useful results as 
signal amplification without the need to do complete simulations. 

Similar amplification may be observed in coupled mechanical 
oscillators. However, here the observed amplification in the energy 
transfer from left (input) to right (output) vortex for (pi,p 2 ,p 3 ) = (1, 
— 1,-1) can be controlled by switching the polarity p 2 from 1 to — 1, 
much like changing the states of a BJT between active and cut-off. 



SCIENTIFIC REPORTS | 4:4108 | DOI: 1 0.1 038/srep041 08 



5 



t = 2.43 ns 




Figure 6 | Fan-out. Stray field distribution in a MVT network. Signal is given only to the left vortex in the middle row. Vortices are marked with their 
respective polarities. A path followed by the antivortex packets taking energy from the left vortex in the middle chain to the left vortex of the top 
chain is marked. The ESDs of left and upper and lower right vortices are given in the inset marked as T, 'Ol' and '02' respectively. 



This can be achieved by using a local magnetic field or a spin-polar- 
ized current; thus, making it a more suitable candidate for integration 
with current electronic technological ecosystem. Moreover the 
observed gain, while remaining constant at B active for low signal 
strength h 0 , drops dramatically for h 0 > 0.25 mT. Thus the output 
will not increase over a certain upper limit (76.61 dB for the MVT 
described here). This is similar to the saturation state of an electronic 
transistor. Direct parallels to all three operational states of a BJT, 
namely active, cut-off and saturation have thus been demonstrated 
for the discussed MVT. Also, both classic transistor operations of 
signal switching and amplification have been described. The depend- 
ence of gain characteristic, in particular 5 active and B max should be 
further explored with different material and geometrical parameters 
of the MVT and the driving signal frequency (dynamic response). 

Our attempt to demonstrate a fan-out operation uncovered that 
the solitons involved in the dynamics do not split easily. This resulted 
in a higher level of asymmetry between different branches of a sym- 
metrical network. This asymmetry was unlikely if the dynamic stray 
field lacked any topologically stable antivortices and treated both the 
branches evenly. This helped us further to validate the importance of 
antivortex solitons in the energy transfer mechanism. More work will 
be needed to demonstrate a successful fan-out operation by consid- 
ering the cascade of antivortices for different network parameters. 

Pinning, which can occur due to fabrication related issues, is 
known to affect the natural frequency of an isolated vortex 48 . 
However, as the dynamics studied here was forced, the observed 
results are expected to remain unaffected unless the pinning poten- 
tial is high enough to change the trajectory of the vortex core or 
the cascade of solitons. A stronger pinning may sometimes occur 
at the Py-air boundary 49 . We have not considered this type of 
pinning here as it can affect the generation and dynamics of the 
stray field which is responsible for some of the reported observa- 
tions. Thus, a different soft-ferromagnet may have to be used if 
pinning becomes an issue. 

While advancing the cause of nano-electronic devices 50 , we also 
hope that these findings will promote the continued search of new 
and improved transistors 51 . For the type of transistor proposed here 
to become technologically viable and competitive, further research 
towards miniaturization and reduction of energy consumption and 
response delay are highly desirable. 



Methods 

Magnetic vortex dynamics was simulated using the finite difference method based 
Landau-Lifshitz-Gilbert (LLG) ordinary differential equation solver called Object 
Oriented Micro -magnetic Framework (OOMMF). Before the dynamics could be 
observed, a magnetic ground state has to be achieved with required vortex core 
polarity and chirality 1,2,6 ' 52 " 57 . This was accomplished by using a pulse field H t = H 0 

exp( — t' 2 ). Here, pi 0 H 0 = 1 T and normalized time f = (t — to) j (V2o^j . t 0 = 75 ps 

and a is the standard deviation of this Gaussian pulse in time whose full width at half 
maximum is 30 ps. Close to the centre of the circular geometries we apply H z = ±H t l 
10 along the Z axis where the sign controls the core's polarity. If the origin of co- 
ordinates is brought to the centre of the vortex then X- and Y- components of fields, 
H x and H y that would produce the desired chirality are given below 

H x =+H t sm(6); 



-_H t cos(9). 



Here, 6 = tan~ 1 (y/x) and the upper or lower signs were chosen for CCW or CW 
chiralities, respectively. It is to be noted that this pulse signal controlled by H t , dies 
down quickly while the magnetic ground state is obtained by running the simulation 
for 40 ns under a high damping (Gilbert damping constant a = 0.95 is used in the 
LLG equation). We have used saturation magnetization, M s = 0.8 X 10 6 A/m, 
exchange constant, A = 13X10~ 12 J/m and zero magneto-crystalline anisotropy. 
During vortex dynamics simulations we reduce a to a more realistic value of 0.008 for 
Py. Magnetization was observed every 10 ps for about 40 ns during dynamics. The 
cell size used during simulation was 5 nm X 5 nm X 40 nm. To ascertain that the 
results discussed here are independent of cell size, Fig. 2 was reproduced using a cell 
size of 2.5 nm X 2.5 nm X 40 nm. The same has been shared as Supplementary Fig. 
F2. Although the form of the peaks have changed to some extent, the relative gains at 
f = fo remain largely unaffected. 

Before we start to explore the dynamics of magnetic vortices, we first need to obtain 
the natural frequencies associated with a single isolated vortex. A broadband 
excitation signal was given to reveal these frequencies. The signal had only X-com- 
ponent, H s x which contained power up tof cut = 45 GHz and depended upon time t as 
given by: 



Ht=H\ 



0 s m(2nf cut (t-t 0 )) 
2nf cut (t-t 0 ) 



Here, fi 0 H° x = 0.05 T and t 0 = 200 ps. 

Upon obtaining the magnetization data from OOMMF, we chose to analyse the 
results by looking at the time evolution of spatial average of normalized X-component 
of magnetization, (m x ) for each vortex, and its corresponding ESD. Normalization is 
done by dividing the X-component of magnetization, M x by M s ; such that m x = MJ 
M s . The Hanning window is used on (m x ) to reduce spectral leakage. The widowed 
data is then zero padded and Fourier transformed to obtain the required ESD, 
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(fh x ) 258 . This is plotted in figures on decibel scale as w H x 20log 10 |ra x |, where a 
window scaling factor of w H = 2 is used for the Hanning window. These ESDs were 
calculated after running the dynamics for over 40 ns, so that any transient vortex core 
dynamics are suppressed and steady state dynamic solutions appear to be more 
prominent in the spectrum. Power spectral density is considered to be more 
desirable in the case of persistent signals. However, here we run the simulations 
for finite amount of time. Also, natural damping ensures that net power input to 
the system becomes zero before the simulations finish. The stray field is also 
obtained from OOMMF during dynamics. The stray field plots shown in Figs. 3 
and 5 and the supplementary movies were created using MATLAB. The contour 
colouring is based on the sum of squares of X- and Y- components of the stray 
field. 
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